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ABSTRACT 

We report two detections of deuterated molecular hydrogen (HD) in QSO absorption-line systems 
at z > 2. Toward J2123-0500, we find iV(HD) = 13.84 ± 0.2 for a sub-DLA with metallicity ~ O.SZq 
and A^(H2) = 17.64 ± 0.15 at z = 2.0594. Toward FJ0812+32, we find Af(HD) = 15.38 ± 0.3 for a 
solar-metallicity DLA with A^(H2) = 19.88 ± 0.2 at z = 2.6265. These systems have ratios of HD to 
H2 above that observed in dense clouds within the Milky Way disk and apparently consistent with a 
simple conversion from the cosmological ratio of D/H. These ratios are not readily explained by any 
available model of HD chemistry and there are no obvious trends with metallicity or molecular content. 
Taken together, these two systems and the two published z > 2 HD-bearing DLAs indicate that HD is 
either less effectively dissociated or more efficiently produced in high-redshift interstellar gas, even at 
low molecular fraction and/or solar metallicity. It is puzzling that such diverse systems should show 
such consistent HD/H2 ratios. Without clear knowledge of all the aspects of HD chemistry that may 
help determine the ratio HD/H2, we conclude that these systems are potentially more revealing of 
gas chemistry than of D/H itself and that it is premature to use such systems to constrain D/H at 
high-redshift. 
Subject headings: ISM: molecules — quasars: absorption lines 



1. INTRODUCTION 

Absorption-line systems in the spectra of high-redshift 
sources such as QSOs and GRB afterglows provide op- 
portunities to examine the interstellar gas in galaxies at 
early cosmic times with samples that are less biased than 
emission-selected samples. These techniques have been 
used to probe the properties of gas clouds at high red- 
shift, where conditions can differ significantly from those 
in the local Universe. Detections of cold molecular gas at 
high redshift are less common than might be expected; 
only about one quarter of DLAs show detectable molec- 
ular hydrogen (Noterdaeme et al. 2008a). Even in GRB- 
DLAs, where the sightline pierces deep inside the ISM 
of the host galaxy, molecular gas does not appear often 
(Tumlinson et al. 2007), though it has been detected in 
one system via H2 and CO absorption (Prochaska et al. 
2009). The opportunity to detect deuterated hydrogen 
molecules (HD) exists in many of these systems where 
H2 is seen, and would in principle provide another use- 
ful constraint on physical conditions. HD is fragile, eas- 
ily photo-dissociated, and is present in the Milky Way 
interstellar medium (ISM) only along reddened sight- 
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lines where the absorbing cloud is dense, dusty, and self- 
shielded. HD may also be an important gas coolant at 
low density in low-metallicity galaxies, and so it is im- 
portant to understand the conditions in which it becomes 
abundant. In this Letter we report two new detections 
of HD in high-redshift QSO absorption line systems: the 
sub-DLA at z = 2.09534 toward J2123-0500 and the 
DLA at z = 2.626 toward FJ0812-H32. These systems 
show H2 and HD absorption but are otherwise quite dif- 
ferent and revealing. 

2. DATA ON J2123-0500 

Data on J2123-0500 were obtained with Keck/HIRES 
in a series of exposures spanning the nights of 17-20 Aug 
2006. A complete analysis of the spectrum is presented 
by Milutinovic et al. (submitted). The effective spectral 
resolution is i? ~ 100, 000 in the reduced spectrum, with 
S/N ~ 15 — 20 pix"^ across the band where we search 
for molecular absorption. The sub-DLA at z = 2.0594 
exhibits logiV(H I) = 19.18 ± 0.15. These data have 
also been used to constrain the temporal evolution of 
the proton-to-electron mass ratio by Malec et al. (2010, 
hereafter MIO). 

We identified the HD lines shown in Figure [1] using the 
laboratory wavelengths of Ivanov et al. (2008) and the 
transition strengths of Abgrall & Roueff (2006). Equiv- 
alent widths W\ (Table 1) were measured by integration 
over the line profile with respect to an iteratively fit- 
ted polynomial continuum. The detected lines were then 
fitted to a single-component curve of growth (COG) pa- 
rameterized by HD column density Af(HD) and the ef- 
fective doppler width bno', the best fit was determined 
to be logiV(HD) = 13.84 and h = 1.2 km s'^. Because 
the detected HD lines lie on the linear and knee por- 
tions of the COG, bHD is poorly constrained at the high 
end; conversely, the column density is well constrained 



TABLE 1 
Summary of two HD systems 



Line 


Ao [A]i 


f' 


W^ [mA]3 
J2123-0050 


W^ [mA]3 
FJ0812+32B 


LO-OR(O) 


1105.84 


0.000743 




7.8 ± 2.0 


Ll-OR(O) 


1092.00 


0.00297 






L2-0R(0) 


1078.83 


0.00675 


4.4 ± 2.8 




L3-0R(0) 


1066.28 


0.0115 


5.6 ± 0.6 


13.1 ± 2.2 


L4-0R(0) 


1054.29 


0.0164 




13.4 ± 1.0 


L5-0R{0) 


1042.85 


0.0206 


8.0 ± 0.6 


12.7 ± 2.2 


L6-0R{0) 


1031.91 


0.0236 






L7-0R{0) 


1021.46 


0.0254 


8.0 ± 1.2 




L8-0R{0) 


1011.46 


0.0262 


8.4 ± 1.3 


13.8 ± 2.5 


WO-OR{0) 


1007.29 


0.0325 


9.7 ± 2.6 


18.4 ± 2.9 


Summary 


logAf(HI) 






19.18±0.15 


21.35±0.10 


logJV(H2) 






17.64±0.15 


19.88±0.10 


logAf(HD) 






13.84±0.20 


15.38±0.20 


bHD [km s-i] 






^9+1.3 


0.85±0.20 



1 Ivanov et al. (2008) 

2 Abgrall & Roueff (2006) 



^ Rest-frame 

to N(HD) > 13.65. As part of their global line-profile 
fitting, MIO obtained logA^(HD) ^ 13.77 and b = 1.9 
km s^^, a consistent result (Figure [T]). 

By contrast with the HD, there is considerable uncer- 
tainty in the H2 column density. MIO assigned two ve- 
locity components to the HD-bearing absorber (1 and 2 
in their Table 3); one of which is narrow and saturated 
in H2. Their evidence for two components was based 
on careful treatments of the profile fits; ultimately two 
overlapping components are preferred by the data once 
continuum and zeropoint adjustments were taken into 
account (see their Section 3 and Figure 5). The column 
densities that minimize x^ for this component are N(0) 
= 15.80, N(l) = 17.52, for a total log iV(H2) = 17.55 and 
bH2 — 1-9 km s^^. These values yield an unphysical rota- 
tional temperature Tqi = 171K/ln(9A^(0)/iV(l)) = -96 
K. Since these lines are highly saturated, much larger 
N(0) is possible with minimal change to the quality of 
the fit. Milutinovic et al. (2010) used only a single com- 
ponent for these profiles in their analysis of the metal 
lines in this system. 

To explore how much higher N(0) could be in this sat- 
urated component, we searched a column-density space 
parameterized by N(0) and N(l). These results appear in 
Figured where we plot contours of x^ for the global fit. 
For this test all other parameters in the global fit were al- 
lowed to settle at their optimal values. The point marked 
"1" is the fiducial fit adopted by MIO, with N(0) = 15.56 
and N(l) = 17.55. We find that the x^ contours extend 
up to N(0) ~ 17 at 99% confidence, well into the range 
where Tqi reaches reasonable physical values of 100 - 200 
K. We therefore adopt a fiducial model that has N(l) = 
17.55 from the MIO profile fit, and N(0) = 16.86 ± 0.24, 
as shown with the range bar in Figure [H This latter 
value splits the difference between the 99% confidence 
interval on N(0) and the line of unphysical Tqij ^nd as- 
signs an error that ranges to these extremes (point "2"). 
This point is within the 99% confidence interval for the 
profile fit and also gives a physically sensible Tqi = 281 
K (with a lower bound of 147 K). This temperature cor- 
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Fig. 1. — H2 and HD lines in the system at 2 = 2.059345 in 
the J2123-0050 sightline. The top panel shows H2 lines from the 
Lyman 2-0 band with the best profile fit overlaid. The middle 
panels show six R(0) transitions of the HD Lyman- Werner bands. 
The lower panel shows the HD parameter space, with the best fits 
from the COG (filled square) and the MIO profile fit (filled circle). 



responds to a thermal bH2 = 1.5 km s^'^ (consistent with 
the profile-fitted b = 1.9 km s~^) and bno = 1-3 km s~^ 
(consistent with the COG-derived value of 1.2). This 
model is consistent with the available line strengths and 
widths for both H2 and HD and becomes our fiducial 
model. The total log7V(H2) = 17.64 ± 0.15 and we find 
logiV(HD) / 2N(H2) = -4.1 ± 0.2. 

We attempted a model of this absorber with a single 
velocity component in the H2 using a COG; for this case 
N(0) = 15.70, N(l) = 16.20, and logiV(H2) = 16.41. 
This fit has the advantage of not requiring any assump- 
tion about the J = column density, but has the dis- 
advantage of being an inferior fit to the data, which 
strongly prefer a two-component fit (MIO). This COG 
has doppler bH2 — 4.6 km s^^, which does not agree as 
well with the buu = 1-9 km s~^ as the profile fit above. 
Most significantly, this much lower N(H2) yields a highly 
anomalous ratio of HD / 2H2 ^ 10^"^. We report this 
"low H2" model here and include this value in later anal- 
ysis for completeness. 
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Fig. 2. — Contours of x^ for the saturated component of H2 
toward J2123-0500. Contours mark the 68, 95, and 99% confidence 
intervals on total N(0) and N(l). Diagonal dotted lines mark the 
locations of Tqi = 100, 200, and 450 K, from right to left. The 
dashed line marks the region where Tqi becomes undefined. Values 
left of this line are excluded on physical grounds. The two adopted 
column densities are marked with circled numerals. 



We do not detect J = 1 HD in this system, but 
we can set a limit. The 5-0R(l) line at 1043.295 A 
falls into a region of clean continuum with good S/N. 
We obtain W\ < 2 mA (3cr), which translates to N(l) 
< 1.5 X lO^'^ cm^^. From this and the measured N(0), 
we infer Tqi < 46 K for this absorber. This temperature 
is lower than the (problematic) Tqi — 292 K for the H2. 
In interstellar H2, Tqi is regarded as a reliable measure 
of the gas temperature only above a density uh ~ 100 
cm~^ at which collisions effectively populate J = 1 from 
J = via proton interchange reactions (Dalgarno, Black, 
& Weisheit 1973). Since HD has a dipole moment, radia- 
tive transitions between J = 1 and J = are permitted 
(as they are not for symmetric H2) and HD can take on a 
lower excitation temperature than H2 even at high den- 
sity. Thus this measurement is probably not a reliable 
indicator of temperature in the HD absorbing cloud. 

Of the Lyman series, only Lya and Ly/3 are covered by 
our data at > 989 A in the rest frame for z = 2.09534. 
The Lya and Ly/3 profiles are damped, and the expected 
D I profile at —82 km s~^ cannot be measured. 

3. DATA ON FJ0812-I-32B 

The data on FJ0812-I-32 were obtained with 
Keck/HIRES in a series of exposures totaling 14,000 s 
over 16-17 Mar 2005. Another exposure was obtained 
on 12 Jan 2008 with the same setup but bluer wave- 
length coverage. The data were reduced using standard 
XIDL reduction packages. The effective spectral reso- 
lution is -R ^ 50, 000 in the reduced spectrum, with 
S/N ^ 40 pix"-'^ across the band where we search for 
H2 and HD, declining to ^ 20 pix^^ in the blue. These 
high-quality data have previously revealed a metal-strong 
DLA (Prochaska, Howk, & Wolfe 2003) with an ex- 
tremely cold gas cloud Jorgenson et al. (2009, hereafter 
J09). 

This system shows six detected HD R(0) lines; the 
other lines missing from Table 1 are blended with other 
features. These lines fall onto the flat portion of a single- 
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Fig. 3.— H2 and HD lines in the DLA at z = 2.626 toward 
FJ0812-I-32. The top panel shows H2 lines from the Lyman 2-0 
band with the best profile fit overlaid. The middle panels show five 
R(0) transitions and one P(l) transition in the HD Lyman- Werner 
bands. The lower panel shows the HD parameter space and the 
best fit from the COG (filled star). The profile-fitted values are 
marked with the cross. 



component COG (Figure [3]), where their column density 
is only moderately constrained. However, the effective 
dopper bHD parameter is well-constrained; the best fit 
is logA^(HD) = 15.38 ± 0.25 and bno = 0.85 ± 0.20 
km s-i from the COG. A profile fit using VPFITEI ob- 
tained log7V(HD) — 15.50 and bno = 1-0, as shown in 
Figure 3. The column densities of H2 in this system were 
also obtained as part of this global profile fit, yielding 
N(J) = 19.81, 19.13, 16.20, 14.86, and 13.97 for J = 0, 
1, 2, 3, 4, respectively, for log N(H2) = 19.88 ± 0.1. The 
best-fit doppler b was held at 1.0 km s~^ for J = and 1 
under the assumption that the bulk of the column den- 
sity in H2 coincides with the cold component of carbon 
and so is thermally broadened at 115 K. These linewidths 
are consistent with the COG and profile-fitted values for 
HD at this temperature. We therefore conclude that the 
HD in this cloud plausibly arises in the same cold gas 
component as the cold C and H2 (J09). 

* http://www.ast.cam.ac.uk/~rfc/vpfit.html 
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Fig. 4. — Comparison of the our two new HD abundances (filled 
circle for J2123-0500 and filled star for FJ0812+32B) relative to H2 
in the Galactic disk from Lacour et al. (2005) (open squares), and 
two DLA systems from Varshalovich et al. (2001, filled triangle), 
and Noterdaeme et al. (2008b, filled diamond). The cross-hatched 
area refers to the axis at right and marks the CMB-derived D/H 
ratio from Dunkley et al. (2009). 

There is also tentative evidence for HD J = 1 in this 
system. Three J = 1 hnes appear as excess absorp- 
tion in blended profiles of unrelated absorption; these 
are Werner O-OR(l) and Lyman 5-0R(l) and 4-0P(l). 
Only the Lyman 4-0P(l) line at 1055.56 A is isolated 
and is marginally detected (see Figure 3). Including all 
four lines in the global fit yields N(l) = 13.50 ± 0.12 for 
bno = 1.0 km s"'^ fixed at the same value as for HD 
J = 0. If we use this tentative value for N(l) we obtain 
Tqi = 26 K, colder than the corresponding C I and H2, 
for which J09 inferred T < 74 K. As for J2123-0500, we 
attribute this low excitation temperature to efficient ra- 
diative decay between J = 1 and J — rather than the 
gas kinetic temperature. 

4. COMPARISONS AND INTERPRETATIONS 

We are now ready to compare these detections to other 
known HD-bearing gas. HD is known to be a component 
of some dense Galactic clouds from detections in sight- 
lines to bright stars in the disk. This population is repre- 
sented in Figure 4 by the sample of Lacour et al. (2005), 
who reported HD toward 17 Galactic OB stars (10 from 
Copernicus, 7 from FUSE). These stars have E(B-V) = 
0.1 — 0.7, with molecular fractions fH2 = 0.06 — 0.7. 

HD has been reported previously in two DLA sys- 
tems. Varshalovich et al. (2001) and, later, Ivanchik 



et al. (2010) reported a detection of HD with N(HD) 
= 3Atli X 10^5 cm-2 and N(H2) = 4.8± 1.0 x lO^^ cm-^ 
in the DLA at z == 2.3377 toward Q1232-H082, which has 
A^(HI) = 20.90 ± 0.08. Also, Noterdaeme et al. (2008b) 
have detected HD in a z = 2.418 DLA toward SDSS 
J1439-H1117 with log N(H2) = 19.38 ± 0.10, log A^(HD) 
= 14.87 ±0.025, and iV(HI) = 20.1 ± 0.1. This HD to 
H2 ratio is also well above the Galactic disk ratio. 

We have plotted two symbols in Figure |4] for the DLA 
toward J2123-0500: the open symbol uses the lower 
COG-fitted measurement for N(H2), while the filled cir- 
cle uses the higher profile fitted value (§2). The error 
bars on the points express the estimated statistical noise 
on each measurement; the line connecting them shows 
the larger systematic uncertainty arising from saturation 
of the H2 lines. With the higher H2 column density, the 
HD to H2 ratio in this system more closely resembles the 
other systems and the cosmological ratio of D/H. 

Our findings in the J2123-0500 and FJ0812+32B sys- 
tems agree in their HD/H2 ratio with the other detec- 
tions, with the striking result that all four high-redshift 
systems lie both well above the Galactic disk value and 
apparently consistent with the cosmological D/H ratio. 
This finding is somewhat puzzling, in light of the diver- 
sity in the other properties of these systems. In their 
analysis of Q1232-I-082, Ivanchik et al. (2010) explain 
this effect by low astration at early times. However, these 
systems show metallicities proxied from S or Zn that are 
between 0.04 Zq and Zq (as shown in the legend of Fig- 
ure |4]). It is puzzling why the infall of gas undepleted in 
D would be arranged just so that the primordial D/H ra- 
tio is maintained in four such diverse systems. It is not 
clear why low astration would be the dominant factor 
in environments that have already achieved a metallic- 
ity comparable to the Milky Way disk, as it is for three 
of the systems. Indeed, the consistency of the HD/H2 
ratio across the factor of 20 in metallicity spanned by 
these systems indicates that metallicity is perhaps only 
a minor factor in determining the abundance of HD. 

Even with four systems that lie along a narrow band 
of HD/H2 ratio well above the Galactic ratio, we can- 
not be sure that this is the general pattern of HD/H2 in 
high-redshift galaxies. HD may have gone undetected in 
other DLAs where the Galactic ratio would place the HD 
below detection limits. Our column density sensitivity 
limits here are N(HD) ~ 10^'^ cm"'^, which gives a limit 
HD/2H2 = 1.1 X 10-5 for J2123-0500 and 6.6 x 10"^ for 
FJ0812-I-32B, well below the Galactic points in the latter 
case but above them in the former. Though in most pub- 
lished cases the data quality is adequate to detect HD in 
the ratios seen here, a systematic survey of H2-bearing 
DLAs might reveal that in some cases a Galactic ratio is 
permitted given the individual detection limits. In light 
of this "detection bias" we cannot yet conclude that the 
anomalous ratios of the four systems known so far are 
representative of the general pattern at high redshift. 

Another striking feature of these systems is that their 
typical molecular fraction is fH2 ~ 0.06 — 0.3. The typi- 
cal value for the Galactic disk clouds is up to ten times 
higher, with most scattering in between fH2 = 0.2 and 
0.6 (Figure[4]). In the Galactic disk, the HD is believed to 
occupy only the dense, self-shielded inner regions of the 
cloud and to be dissociated easily in the low-extinction 



photodissocation regions surrounding the cloud (Lacour 
et al. 2005). By this hne of reasoning, clouds with lower 
molecular fractions should have more of their HD dis- 
sociated and should lie below the Galactic disk points 
rather than above. This effect is particularly notable for 
the J2123-0500 system, which does not qualify as a DLA 
and has the lowest fH2 — 0.028. The mystery is then 
deepened, in the sense that we must explain a higher HD 
to H2 ratio than in the Milky Way in clouds of similar 
metallicity but lower total molecular abundance. 

Dust is another factor that may help promote the for- 
mation of HD in these systems. All four of these HD- 
bearing systems show evidence of dust depletion based 
on dust-sensitive abundance ratios. The system toward 
J2 123-0500 shows ionization-corrected abundance ratios 
Si and Fe depleted by approximately 0.5 dex with re- 
spect to S and N (see table 4 of Milutinovic et al.). The 
F0812-^32B system shows solar [S/Si] but [Zn/Fe] and 
[Cr/Zn] that indicates a high depletion of Fe (J09). The 
2175 A dust feature is not detected in either spectrum. 
Unfortunately, beyond these abunance ratios no more 
quantitative information has been obtained on the exact 
dust-to-gas ratios or dust properties of these two systems, 
so it is not known whether the composition and size dis- 
tribution of the dust is similar to that of the Galactic 
disks. We therefore cannot correlate the HD properties 
of these systems against dust in any more detail. 

Using a master-equation approach, Lipshtat, Biham, 
& Herbst (2004) studied the relative formation rates of 
H2 and HD under the assumption that D atoms stick to 
dust grain surfaces for slightly longer than H. They have 
found that for a certain range of conditions, the total 
formation rate of HD can rise from ~ 10~^ to > 1Q~^ 
times that of H2 as the dust temperature rises from < 15 
K to 20 K. The effect is strongest for small grains with 
a small count of binding sites to which H and D attach. 

Since this model describes the formation rates but not 
the relative photo-dissociation rates (which are greatly 
influenced by self-shielding for H2 and much less so for 
HD), and since it possesses at least three free parame- 
ters for density, grain size, and grain temperature, we 



cannot directly use our data to constrain these models. 
We can speculate that small grains, the right range of 
dust temperatures, or both, may together establish the 
observed molecular abundance ratios seen in these four 
z > 2 systems. However, with the relative rates of H2 
and HD production ranging over a factor of 100 between 
Tdust < 15 K and 20 K for a range of grain sizes, it is 
not clear why the four systems found so far have scatter 
so closely around the cosmological D/H. 

We also compare our HD results against the D/H ra- 
tios seen in high-redshift DLAs as compiled by Pettini 
et al. (2008), represented by open diamonds in Figure ID 
We note that the D/H measurements are quite precise 
and show a low scatter consistent with the CMB-derived 
value. In light of the precision and consistency of these 
measurements, it seems unnecessary to rely on HD/H2 
as a measure of cosmological D/H when the chemical 
influences on HD are so uncertain. 

In light of the diversity of these systems in metallicity 
and molecular fraction, and with the discussion above 
in mind, we conclude that the balance of chemical re- 
actions, dust properties, and local conditions are highly 
uncertain in these systems. Following this, it seems that 
these systems pose a puzzle about interstellar chemistry 
more effectively than they solve one about D/H. If prop- 
erly understood these systems are likely to be far more 
revealing of these uncertain chemical and environmental 
effects than they are of the primordial abundances, which 
in any case appear well constrained from direct measure- 
ments of D/H. We can hope the discovery of more such 
systems and additional theoretical work will be able to 
disentangle these effects and resolve these doubts. 
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